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Abstract 

We present the calculation of scalar and pseudoscalar Higgs production in 
association with a top-antitop pair to the next-to-leading order (NLO) accu- 
racy in QCD, interfaced with parton showers according to the MC@NLO 
formalism. We apply our results to the cases of light and very light Higgs 
boson production at the LHC, giving results for total rates as well as for 
sample differential distributions, relevant to the Higgs, to the top quarks, 
and to their decay products. This work constitutes the first phenomeno- 
logical application of aMC@NLO, a fully automated approach to complete 
event generation at NLO in QCD. 
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1. Introduction 



Establishing evidence for the Higgs boson(s), i.e., the scalar remnant(s) 
of the Englert-Brout-Higgs mechanism [l|, 0, 0] in the standard model and in 
extensions thereof, is among the most challenging goals of the LHC experi- 
ments. A coordinated theoretical/experimental effort in the last years has 
led to a number of remarkable achievements in the accuracy and usefulness 
of the available theoretical predictions, and in the role these play in current 
analysis techniques jij. 

Depending on mass and couplings, Higgs bosons are produced and even- 
tually decay in a plethora of different ways, leading to a wide range of 
signatures. In most cases, signals are difficult to identify because of the 
presence of large backgrounds, and reliable predictions are necessary firstly 
to design efficient search strategies, and secondly to perform the correspond- 
ing analyses. A particularly challenging scenario at the LHC is that of a 
standard-model light Higgs, m H < 130 GeV. In this case, the dominant de- 
cay mode is into a bb pair, which is however completely overwhelmed by the 
irreducible QCD background. A possible solution is that of considering the 
Higgs in association with other easier-to-tag particles. An interesting case 
is that of a top-antitop pair, since the large Yukawa coupling ttH , and the 
presence of top quarks, can be exploited to extract the signal from its QCD 
multi-jet backgrounds. Unfortunately, this production mechanism is also 
plagued by large backgrounds that involve a ti pair, and hampered by its 
rather small rates, and thus turns out to be difficult to single out. Several 
search strategies have been proposed, based on different decay modes: from 
bb which leads to largest number of expected events, to the more rare but 
potentially cleaner rr [fj , H / W / **- ) and 77 0] final states. All of them are 
in fact very challenging, and dedicated efforts need be made. For example, 
recently it has been argued that in the kinematical regions where the Higgs 
is at quite high transverse momentum the bb pair would be merged into one 
"fat" jet, whose typical structure could help in discriminating it from QCD 
backgrounds @, 01 (boosted Higgs scenario). 

It is then clear that accurate and flexible simulations, for both signals 
and backgrounds, can give a significant contribution to the success of any 
given analysis. Predictions accurate to NLO in QCD and at the parton 
level for ttH hadroproduction have been known for some time 0, Q 0, 



13l . Il4j . Il5| , and recently confirmed by other groups [16|, ll7j • As for the most 



relevant background processes to the Higgs decay mode into bb, NLO cal- 



culations for ttbb 18 



and ttjj [2l| are available in the literature. In 
2 



this work, we extend the results for the signal to computing the associated 
production ttA of a pseudo-scalar Higgs boson. All aspects of the calcu- 
lations we present here are fully automated. One-loop contributions have 
been evaluated with MadLoop [17], that uses the OPP integrand reduction 



method |22J as implemented in CutTools [23|. The other matrix-element 



contributions to the cross sections, their phase-space subtractions accord- 



ing to the FKS formalism [241 ] . their combinations with the one-loop results 



and their integration are performed by MadFKS [25|. The validation of 



MadLoop and MadFKS in the context of hadronic collisions has been 



presented in Ref. [17J. For the sake of the present work, we have also per- 
formed a dedicated comparison with the results of Ref. [4j for the total tiH 
cross section, and found agreement at the permille level for several Higgs 
masses. 

We have also matched our NLO results with parton showers using the 
MCONLO method 0. This matching procedure has also been completely 
automated, and this work represents the first application of the MC@NLO 
technique to non-trivial processes which were previously available only at 
fixed order and at the parton level - in other words, to processes not already 
matched to showers by means of a dedicated, final-state-specific, software. 
What said above also implies that our results are the first example of NLO 
computations matched to showers in which all ingredients of the calculation 
are automated, and integrated in a unique software framework. 

We remind the reader that the structure of the MC@NLO short-distance 
cross sections is the same as that of the underlying NLO computation, ex- 
cept for a pair of extra contributions, called MC subtraction terms. These 
terms have a factorised form, namely, they are essentially equal to the 
Born matrix elements, times a kernel whose main property is that of being 
process-independent. This is what renders it possible the automation of the 
construction of the MC subtraction terms, and ultimately the implemen- 
tation of the MC@NLO prescription. We call aMC@NLO the code that 
automates the MC@NLO matching, and we defer its detailed presentation 
to a forthcoming paper 27| . aMC@NLO uses MadFKS for phase-space 
generation and for the computation of the pure-NLO short distance cross 
section of non-virtual origin, and on top of that it computes the MC sub- 
traction terms. One-loop contributions may be taken from any program 
which evaluates virtual corrections and is compatible with the Binoth-Les 
Houches format [28]; as was said before, we use MadLoop for the predic- 
tions given in this work. The resulting MC@NLO partonic cross sections 
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are integrated and unweighted by MINT |2j|, or by BASES/SPRING [3C 
aMC@NLO finally writes a Les Houches file with MC-readable hard events 
(which thus includes information on particles identities and their colour con- 
nections). 

2. Results at the LHC 

We present selected results for total cross sections and distributions rel- 
evant to ttH/ttA production at the LHC in three scenarios: 

I. Scalar H, with tuh = 120 GeV; 

II. Pseudoscalar A, with uia = 120 GeV; 

III. Pseudoscalar A, with = 40 GeV; 

where the Yukawa coupling to the top is always assumed SM-like, ytj \/2 = 
m t /v. 

The three scenarios above allow one to compare the effects due the dif- 
ferent parity of the Higgs couplings on total rates as well as on differential 
distributions. In this respect, it is particularly interesting to consider the 
situation in which the Higgs boson is light and pseudoscalar, as is predicted 



in several beyond-the-standard-model theories (see e.g. Refs. |3ll. l32l. l33|). 
The main purpose of this section is that of studying the impact of QCD 
NLO corrections at both the parton level and after shower and hadronisa- 

tion. For the nu m erica. analysis we choose „, = „.= (< m <<'f , 

where rriT = \/ m 2 + p% and mF t ole = mf IS = 172.5 GeV. We have used LO 
and NLO MSTW2008 parton distribution functions for the corresponding 
cross sections. The parton shower in aMC@NLO has been performed with 
Fortran Herwig Q@,0, version 6.520 0. 



The predicted production rates at the LHC running at yfs = 7 and 14 
TeV are given in Table [1] where, for ease of reading, we also show the fully 
inclusive K-factor. As far as differential distributions are concerned, we 



1 These integrators have been modified by us, in order to give them the possibility of 
dealing with both positive- and negative-weighted events. 

2 We remind the reader that the MC@NLO formalism has been employed to match 
NLO results with Herwig++ 37 1 and, to a lesser extent, with Pythia HI] (see Ref. [39[ 
and Ref. respectively). The automation of the matching to these event generators is 
currently under way. 
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restrict ourselves to the 7 TeV LHC, and begin by studying a few fully- 
inclusive ones (see Figs. dM]). We then consider a "boosted" case, i.e. apply 
a hard cut on the transverse momentum of the Higgs (see Figs. Eland [HI). 
Finally, in Figs. [7] and |8] we present our aMC@NLO predictions for corre- 
lations constructed with final-state B hadrons, which may or may not arise 
from the decays of the Higgs and/or of the tops (see a discussion on this 
point later). 

We first note a very interesting feature of Fig. [TJ the px distributions 
corresponding to the three different scenarios, while significantly different at 
small transverse momenta, become quite close to each other at higher values. 
This is expected from the known pattern of the Higgs radiation off top 
quarks at high in both the scalar and the pseudoscalar cases 
This difference is not affected by NLO corrections, and could therefore be 
exploited to identify the parity of the coupling at low pr- On the other 
hand, the independence of the parity and masses of the pr distributions at 
high values implies that the boosted analyses can equally well be used for 
pseudoscalar states. 

In general, we find that differences between LO and aMC@L(j|f|, and 
between NLO and aMC@NLO, are quite small for observables involv- 
ing single-inclusive distributions, see Figs. [T]|3j The same remark applies 
to the comparison between LO and NLO, and between aMC@LO and 
aMC@NLO. However, if the cut p" /A > 200 GeV is imposed (boosted 
Higgs analysis), differences between LO and NLO (with or without showers) 
are more significant, and cannot be approximated by a constant if-factor. 

As is obvious, the impact of the shower is clearly visible in the three- 
particle pT{tiH/tiA) distribution of Fig. HI This observable is infrared- 
sensitive at the pure-NLO level for px — > 0, where it diverges logarithmi- 
cally. On the other hand, the predictions obtained after interfacing with 
shower do display the usual Sudakov suppression in the small-pr region. 
At large transverse momenta the aMC@NLO and NLO predictions coin- 
cide in shape and absolute normalisation, as prescribed by the MC@NLO 
formalism. 



3 We call aMC@LO the analogue of aMC@NLO, in which the short-distance cross 
sections are computed at the LO rather than at the NLO. Its results are therefore equiv- 
alent to those one would obtain by using, e.g., MadGraph/MadEvent [43} interfaced 
to showers. 
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Cross section (fb) 


Scenario 




7 TeV 




14 TeV 




LO 


NLO 


.fT-factor 


LO 


NLO 


.fT-factor 


I 


104.5 


103.4 


0.99 


642 


708 


1.10 


II 


27.6 


31.9 


1.16 


244 


289 


1.18 


III 


69.6 


77.3 


1.11 


516 


599 


1.16 



Table 1: Total cross sections for ttH and ttA production at the LHC (y/s — 7, 14 TeV), 
to LO and NLO accuracy. The integration uncertainty is always well below 1%. Scale 
choices and parameters are given in the text. 

In our Monte Carlo simulations, we have included the t — > e + ub, t — > 
e~~vb, and H — >■ bb decays at LO with their branching ratios set to one9. 
After showering, the b quarks emerging from the decays of the primary par- 
ticles will result into ^-flavoured hadrons. As prescribed by the MC@NLO 
formalism, the showering and hadronisation steps are performed by the 
event generator the NLO computation is matched to, i.e. Herwig in this 
Letter. The parameters that control hadron formation through cluster de- 
cays are set to their default values in Herwig [36j. Additional 6-flavoured 
hadrons may be produced as a consequence of g — > bb branchings in the 
shower phase. For example, for scalar Higgs production at 7 TeV, about 
2.7% and 0.5% of events have six and eight lowest-lying B hadrons respec- 
tively. In our analysis, we have searched the final state for all lowest-lying 
B hadrons, and defined two pairs out of them, a) The pair with the largest 
and next-to-largest transverse momenta; b) the pair with the largest and 
next-to-largest transverse momenta among those B hadrons whose parent 
parton was one of the b quarks emerging from the decay of the Higgs (there 
are about 0.2% of events with four or six B hadrons connected with the 
Higgs). The definition of b) relies on MC truth (and in all cases we assume 



4 We have neglected production angular correlations [44[ , as these are expected to have 
a minor impact for the kind of processes and observables we consider here. As usual when 
matching fixed order calculations to parton showers, colour information is transferred in 
the large- Nc limit. The 6-quark mass in the top and Higgs decay products has been set 
to the Herwig default, 4.95 GeV. 
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100% tagging efficiency), but this is sufficient to study the basic features of 
final-state B hadrons. 

In Figs. [7] and [S] we plot the pair invariant mass (tubb) an d the r] — ip dis- 
tance (ARbb) correlations between the I?-hadron pairs defined as explained 
above. The effects of the NLO corrections to tiH/tiA are, in general, mod- 
erate. A cut of 200 GeV on the pt of the Higgs is seen to help discriminate 
the B hadrons arising from the Higgs from those coming either from top 
decays, or from the shower. The shapes of the distributions are similar be- 
tween scenarios I and II while, due to the lower Higgs mass, the tubb and 
ARbb histograms peak at lower values in the case of a pseudoscalar A with 
m A = 40 GeV. 




100 200 300 400 500 600 



Pt HA [GeV] 

Figure 1: Higgs transverse momentum distributions in tiH/tiA events at the LHC (y/s=7 
TeV), with aMC@NLO in the three scenarios described in the text: Scalar (blue) and 
pseudoscalar (magenta) Higgs with ra H / A — 120 GeV and pseudoscalar (green) with 
m A = 40 GeV. In the lower panels, the ratios of aMC@NLO over LO (dashed), NLO 
(solid), and aMC@LO (crosses) are shown for each scenario. 
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Figure 2: Same as in Fig. [TJ for the pr of the top quark. 
3. Conclusions 

Accurate and flexible predictions for Higgs physics will play an impor- 
tant role in understanding the nature of the EWSB sector in the standard 
model and beyond. In this Letter we have presented the results at NLO 
in QCD for (scalar and pseudoscalar) Higgs production in association with 
a top-antitop quark pair, both with and without the matching to parton 
showers. Our approach is fully general and completely automated. A sim- 
ple study performed on key observables involving the Higgs, the top quarks, 
and their decay products shows that while changes in the overall rates can 
be up to almost +20% (for the pseudoscalar states) with respect to LO 
predictions, in general the shapes of distributions are mildly affected for a 
light SM Higgs. Significant changes, however, can be observed in the case 
of a light or very light pseudoscalar state. 

The kernels of MC subtraction terms defined in the M C@NLO formal- 
ism, although process-independent, do depend on the specific event gener- 
ator one adopts for the shower phase. In other words, each event generator 
requires a set of MC subtraction terms, which are computed analytically. 
Results are now available for the cases of Fortran Herwig, Herwig+ + , 
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Figure 3: Same as in Fig. [TJ for the invariant mass of the top-antitop pair. 

and Pythia6; those relevant to the former program have been used to ob- 
tain the predictions presented here, while those relevant to the latter two 
codes are presently being automated and tested against known benchmarks. 

We conclude by pointing out that work is in progress to make the use 
of aMC@NLO for tiH/tiA production and for other processes publicly 
available at http : / / amcatnlo . cern . ch. 
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Figure 7: Invariant mass distributions of the _B-hadron pairs denned as a) (red) and b) 
(blue) in the text. The results obtained by imposing p T > 200 GeV (magenta and 
cyan, respectively) are also displayed. Solid histograms are aMC@NLO, dashed ones 
are aMC@LO. 14 
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